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Abstract. I report on a new method for measuring deviations from ax- 
isymmetry in the velocity fields of gas disks. The method is an extension 
of the method that Franx, van Gorkom and de Zeeuw developed for a 
single orbit gas ring. The measurement is based upon a higher order 
harmonic expansion of the full velocity field. Epicyclic theory is used to 
calculate the effects of a perturbation in the potential on these harmonic 
terms. It is shown how the s\ and S3 harmonics can be used to measure 
€r sm(2(p b s ), where e# is the elongation of the potential and tf> b s is an 
(unknown) viewing angle. The advantage of this method over previous 
attempts to derive the elongation of dark matter halos is that, by using 
HI, we can probe the potential at radii beyond the stellar disk, into the 
regime where the dark matter is thought to be the dominant dynamical 
component. 

As a demonstration, I applied this method to HI velocity fields of 
two spiral galaxies, NGC 2403 and NGC 3198. NGC 3198 shows a very 
small e^j sin(2<^ & s ), which suggests that the galaxy is very nearly axisym- 
metric. NGC 2403 shows a larger €r sin(2<^> b s ), varying in a systematic 
way between and 0.1 and is probably influenced by spiral arms. The 
current data suggests that spiral galaxies are close to axisymmetry, but 
a larger sample is needed to quantify this statement. 



1. Introduction 

It is generally accepted that disk galaxies have massive dark halos, but little 
is known about the shape of the dark matter halos. Dark halos were generally 
modeled as being spherical until Binney (1978) argued that the natural shape 
of dark halos is triaxial. If halos are indeed triaxial, the disk of the galaxy will 
reside in one of the principal planes of the dark halo. In a triaxial halo the axes 
obey the relation a > b > c. I will assume that the disks of spiral galaxies lie 
within the a, b plane. 

Triaxial dark matter halos occur naturally in cosmological N-body simula- 
tions of structure formation in the universe (Katz & Gunn 1991, Dubinski & 
Carlberg 1991, Dubinski 1994), but the exact distribution of the halo shapes is 
as yet uncertain. Especially the inclusion of baryons in the simulations tends 
to make the halos more oblate (Dubinski 1994). So the observed distribution 
of the shapes of dark matter halos will be a powerful constraint on scenarios of 
galaxy formation and subsequent evolution. 



Measuring the shapes of dark halos can be split into two parts: measurement 
of the ratio c/a, i.e. the flattening perpendicular to the plane of the disk, and 
measurement of the intermediate to major axis ratio b/a, i.e. the elongation in 
the plane of the disk. The axis ratio c/a can be measured using, for instance, 
polar ring galaxies (e.g. Whitmore, McElroy & Schweizer 1987, Sackett & Sparke 
1990, Sackett et al 1994), kinematics of halo stars of the Milky Way (van der 
Marel 1991) or the flaring of Hi-disks (Oiling, these proceedings, Sicking 1996). 
All these measurements seem to indicate that c/a < 1, in other words, a spherical 
halo is excluded by these measurements. 

The elongation of the dark halo (b/a) is a similarly uncertain parameter. 
Attempts to measure this quantity include studies of the inclination distributions 
of large samples of spiral galaxies and fitting models with different elongations 
to them (Lambas 1992, Binney & de Vaucouleurs 1981, Fasano et al 1992), 
detailed observations of the old stellar disks of spiral galaxies (Rix & Zaritsky 
1995), kinematics of the milky way (Kuijken & Tremaine 1994) and the scatter 
in the Tully-Fisher relation (Franx & de Zeeuw 1992). All these measurements 
indicate that b/a ~ 0.9. 

The only direct measurement of the elongation of the potential of a disk 
galaxy so far is of the SO galaxy IC 2006 (Franx, van Gorkom & de Zeeuw 1994, 
hereafter FvGdZ), which contains a large HI ring in the plane of the disk. Using 
information on both the kinematics and the geometry of the ring, FvGdZ were 
able to measure the elongation of the potential at the location of the ring. The 
measured elongation was consistent with zero. The method used by FvGdZ 
was based on the assumption that one could use epicycle theory to predict the 
geometry and velocity field of the ring in a mildly perturbed potential. If the 
external potential is elongated, the velocity variation along the ring will not 
be precisely sinusoidal, but higher harmonic terms will be superimposed on it. 
Making a harmonic expansion of the velocity field of the ring and interpreting 
the measured harmonics within the framework of epicycle theory, given the ring 
geometry, the elongation of the potential at the position of the ring can be found. 

I have extended the FvGdZ formalism for measuring the elongation of the 
potential of a single orbit gas ring to the case of a slightly non-axisymmetric 
gas disk, which may contain spiral-like perturbations. This analysis assumes 
a stationary perturbation and closed stable orbits. Therefore, applicability to 
non-linear phenomena like spiral arms is limited. But a small global elongation 
of the overall potential, as is the case with a triaxial dark matter halo, can 
be analysed with this method. The advantage of this method over previous 
attempts to derive the elongation of dark matter halos is that, by using HI, we 
can probe the potential at radii beyond the stellar disk, into the regime where 
the dark matter is thought to be the dominant dynamical component. 

2. Results from perturbation theory 

Suppose the (cold) Hi-disk resides in a potential 

V(R, 0) = V (R) + V m (R) cos(m^ + <p m {R)), (1) 

where R and <f> are polar coordinates in a frame that rotates with the perturba- 
tion, at a pattern speed f2 miP . Here, Vo(R) is the unperturbed potential, V m (R) 
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the amplitude of the perturbation and ip m (R) its phase. The perturbation is 
assumed to be stationary. One can calculate the possible closed loop orbits of 
the gas in this potential (analogous to Binney & Tremaine, 1987, p. 146). Sub- 
sequently one can determine the velocity field generated by these orbits. Now 
assume that this velocity field is observed under viewing angles (4>obs,i), where 
4> b s is the angle between the line where eft = and the projection of the line-of- 
sight onto the plane of the disk and i is the inclination of the disk. If one then 
introduces the azimuthal angle tfj = <p — <j) ohs + tt/2 (the angle in the plane of the 
disk that is zero on the line of nodes), it is possible to show that the velocity 
field projected on the sky (hereafter "line-of-sight velocity field" or "l.o.s. veloc- 
ity field") has the following form (express the line-of-sight velocity as a Fourier 
expansion vi os = J2 n °n cos n V> + s n sin nip): 

vios = v* [ci cos if) + s m -i sin(m - l)ip + c m -\ cos(m - l)ip 

+ s m+ i s'm(m + 1)V> + c m+ i cos(m + l)ip] , (2) 

with v* = v c sini (v c is the circular velocity) and coefficients c x ,s x that depend 
on m, (fiobs, (f(R), Vo(R), V m (R) and TOiP . Now the measurable parameters 
(v*,c x , s x ) are expressed in terms of internal parameters for a potential perturbed 
by a single Fourier term. 

From equation (|2|) we can conclude that if the potential has a perturbation 
of Fourier number m, the l.o.s. velocity field contains m — 1 and m + 1 Fourier 
terms. Qualitatively, this conclusion was also inferred by Canzian (1993). 

In order to measure the harmonic terms in the l.o.s. velocity field, we first 
fit a set of tilted-rings to the l.o.s. velocity field using standard tilted-ring fitting 
routines. Then the velocity along each ring is decomposed into its harmonics. In 
general, one does not a priori know the inclination i, position angle T and centre 
(^0) Ho) °f a g a l ax Y- Since the l.o.s. velocity field itself is used to determine these 
parameters, it is clear that this can influence the resulting harmonic terms. 
FvGdZ showed that this fitting of the ring parameters from the l.o.s velocity 
field will affect the expansion given in equation @ due to a possible difference 
between the best fitting parameters (i and T) and the true parameters. It can be 
shown that only for m = 2 there will be first order differences between the best 
fitting and true internal parameters, whereas for m ^ 2, equation (Q) predicts 
the correct form of the l.o.s. velocity field to second order. In the m = 2 case, 
a misfit of the kinematic centre will result in additional co,S2 and C2 terms in 
the harmonic expansion and misfitting of the inclination or position angle will 
affect the ci,si and 03,53 terms. 

In the case of a globally elongated potential, as caused by for instance a 
triaxial halo, we have m = 2, 0,2,p ~ and <f2{R) ~ const. Note that here we 
are in the m = 2 case, where the best fitting i and T are not equal to the true i 
and r. Assuming a flat rotation curve, one finds that after a tilted-ring fit the 
l.o.s. velocity field has the following form (as derived earlier by FvGdZ): 

vios = ci cos i]} + si sin ip + C3 cos 30 + S3 sin 3tp, (3) 

with 

ci = v c smi(l - ±e v cos2(p obs ), 
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where q = cos 2, ip is the angle in the plane of the ring is zero on the 
apparent major axis of the ring), €r is equal to the elongation of the potential 
and e v = 2e#. We see that from the s± and S3 terms we are able to derive 
e^j s'm(2(fi b s ) , the quantity that I will try to measure. Unfortunately, one cannot 
determine 4>obs separately, so the combination of viewing angle and ellipticity 
is all one can determine for an individual galaxy. Global ellipticity will create 
s\ and S3 terms that are constant with radius (and thus €r sin(20 o ;, s ) will be 
constant as a function of radius as well), whereas spiral- like perturbances will 
cause si and S3 terms that change sign as a function of radius and en sin(20 o & s ) 
will wiggle. Therefore, I will interpret an average value of €r sin(20 o { )S ) that is 
offset from zero as due to ellipticity of the disk. 

Strong warping will invalidate the approach discussed above. 



3. Some example velocity fields 

In order to clarify the equations presented in section 2, we present a number 
of example potential perturbations and corresponding residual velocity fields 
in figure [1]. The harmonic terms of the different potential perturbations range 
from m = 1 to m = 4. For the phase of the fields, ip m {R), we chose a logarithmic 
spiral. The amplitude of the perturbation has been taken constant throughout 
the field. The l.o.s. velocity fields are then created using equation @ (these are 
not shown). A tilted-ring fit was then made to these l.o.s. velocity fields and 
subsequently the harmonic fit was made, revealing the two harmonic components 
that were hidden in them. 

In figure [1] we see from left to right 

1. The m term potential perturbation. 

2. The residual velocity field (velocity field minus fitted circular velocity) 
caused by this potential perturbation (i.e. the m — 1 plus m + 1 terms 
together). 

3. The m + 1 component of the residual field. 

We see that indeed an m term in the density causes m — 1 and m + 1 terms 
in the kinematics of the gas. In these examples, the m — 1 term is dominant. 

Equation (2) predicts that the C3 term should be fitted to zero by the 
tilted-ring fit in the case of a m = 2 perturbation in the potential. Thus, in 
the m + 1 component of the residual field, no spiral-like structure should be 
visible any more. In figure [1], we see that this is indeed the case. Instead of 
spiral structure, we see radially alternating positive and negative contributions, 
caused solely by the S3 (R) term. If the potential would have contained a global 
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Potential perturbation residual field in ] component 




Figure 1. Potential perturbations and velocity fields. We created a 
potential with a perturbation as given in equation (|]) with a single m 
(ranging from 1 to 4, shown in column 1). Subsequently we calculated 
the velocity field corresponding to this potential using equation (2) and 
made a tilted-ring fit to this field. This gave us a residual velocity field 
(column 2). This residual field is dominated by the m — 1 component. 
A harmonic fit to the velocity field revealed that the residual field 
contained the m+ 1 term as well (column 3). 

intrinsic ellipticity, m = 2,(f2(R) = const, the si and S3 terms caused by this 
elongation would have been constant as a function of radius. Therefore, the 
measured S3 would have been constant and non-zero as well and no radially 
alternating pattern for S3 would have been visible in this plot. 

Since in the fitting procedure the centre was kept fixed, we also see a two- 
armed spiral in the m + 1 residual field of the m = 1 potential perturbation. 
If we would have taken the centre as a free parameter in the fit, it would have 
drifted in such a way as to make the C2 and S2 terms disappear. 

4. Two test cases: NGC 2403 and NGC 3198 
4.1. Data description 

NGC 2403 is a nearby (3.25 Mpc; Begeman, 1987) Sc(s)III galaxy (Sandage & 
Tammann, 1981). It has been observed with the WSRT, 4 x 12 hours (Sicking, 
1996). The data have been smoothed to a circular beam of 13 arcseconds. The 
second relatively nearby (9.4 Mpc; Begeman, 1987) spiral galaxy we will examine 
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Figure 2. Harmonic decomposition of the velocity field of NGC 3198. 
Only the first four terms of the expansion are shown, since higher order 
terms have low amplitudes. The c\ term gives the circular velocity, i.e. 
the rotation curve. From the S\ and S3 terms €r(R) sin(20 o f, s ) will be 
measured. The C2, S2 terms are zero if the kinematic centre of a ring is 
left free in the fit. Since the same (fixed) centre is taken for all rings 
here, the 02,^2 terms vary. Especially in the inner parts the C2 term 
is relatively large, indicating that NGC 3198 might be somewhat lop- 
sided. The C3 term is zero if the inclination is fitted correctly (equation 
4), which is obviously the case here. There is also some power in the 
04,54 modes. This may be caused by some three- or five-armed spiral 
structure (see figure 1). Visual inspection of the HI surface density 
map shows a three-fold spiral-like structure. The arrows denote the 
Holmberg radius. 
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Figure 3. cr(R) sin(2c/> & s ) as a function of radius. The effect in 
NGC 2403 may be caused by spiral arms. The average value is used as 
an upper limit to the halo ellipticity. The measurement for NGC 3198 
is relatively straight within errors and very low. The dotted line is 
the average value of €r(R) sin(2</> 6 s ), the arrow denotes the Holmberg 
radius. 
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is the Sc(rs)I-II (Sandage & Tammann, 1981) galaxy NGC 3198. This galaxy 
was observed by Sicking (1996) with the WSRT, also for 4 x 12 hours. The 
observations have a circular beam with FWHM of 18 arcsec. Both galaxies 
are well suited for our type of analysis: they have large extents on the sky 
(Holmberg dimensions: NGC 2403: 29'. x 15'.0 and NGC 3198: ll'.9 x 4'. 9), 
have a favourable inclination (NGC 2403: i m 62°, NGC 3198: % m 71°), the data 
have high signal/noise, the beams are small and both galaxies show no obvious 
warp. The residual maps of the galaxies show only small (typically < 10 km s _1 ) 
residuals, but they appear to contain systemic structures. 

4.2. Measuring the elongation 

The first step in our procedure is to fit concentric tilted-rings to the l.o.s. velocity 
fields of these galaxies. The width of each ring is taken to be about twice the 
beam size, to be sure that all the individual rings are independent. Table Q 
lists the radially averaged ring parameters found by the tilted-ring fit. 



Parameter 


NGC 2403 


NGC 3198 


position angle 
inclination 
systemic velocity 
e R s'm(2(j) obs ) 


124?1±0?09 
61?5 ± 0?20 
133.3 ±0.10 km s" 1 
0.064 ± 0.003 


215?93±0?13 
70?6 ± 0?21 
660.1 ±0.16 km s" 1 
-0.019 ± 0.003 



Table 1. Radially averaged ring parameters for NGC 2403 and 
NGC 3198 



After the tilted-ring fit a harmonic fit was made to the velocity fields along 
each individual ring. Harmonic terms were measured up to ninth order. In 
figure [2], the result for the first four harmonics of NGC 3198 are shown as 
an example. The C3 term is zero everywhere in both galaxies, indicating that 
the inclination is correctly fitted for all rings. Furthermore, we find that in 
both galaxies harmonics c n ,s n with n > 5 are consistent with zero. In the 
lower harmonics, clear systematic trends (as a function of radius) are visible. 
Using si(R) and ss(R), we are able to measure €r(R) sin(20 o & s ). The results of 
these measurements are given in figure [3], where we see that e R (R) sm(2(p b s ) 
strongly wiggles as a function of radius for NGC 2403, whereas for NGC 3198 
it is fairly constant within errors. Furthermore, note that en(R)sm(2cp bs) > 
for NGC 2403 at all radii. This indicates that, although spiral arms are 
present that cause the wiggling, there may also be a global elongation present, 
lifting €r(R) sin(20 o b s ) from an average of to an average of 0.064 ± 0.003. For 
NGC 3198 the average value of e R sm(2(j) obs ) is very low: -0.019 ± 0.003. 

Since the effect of ellipticity and spiral arms cannot be separated unambigu- 
ously, these measurements suggest that spiral galaxies are close to axisymmetric. 
But more galaxies are needed in order to quantify this result (Schoenmakers et 
al, in preparation). 
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